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The nature of the Mott-Hubbard metal-insulator transition in the infinite-dimensional Hubbard model is 
investigated by Quantum Monte Carlo simulations down to temperature T = W/UO {W = bandwidth). 
Calculating with significantly higher precision than in previous work, we show that the hysteresis below 
TiPT — 0.022W, reported in earlier studies, disappears. Hence the transition is found to be continuous 
rather than discontinuous down to at least T = 0.325Tipt. We also study the changes in the density of 
states across the transition, which illustrate that the Fermi liquid breaks down before the gap opens. 
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The explanation of the nature of the Mott-Hubbard 
metal-insulator transition, i.e., the transition between a 
paramagnetic metal and a paramagnetic insulator, is one 
of the classic and fundamental problems in condensed 
matter physics ||l|,g] . Metal-insulator transitions of this 
type are, for example, found in transition metal oxides 
with partially filled bands near the Fermi level. For 
such systems band theory typically predicts metallic be- 
havior. The most famous example is V2O3 doped with 
Cr l^^l^. In particular, in (Vq. 96010.04)203 the metal- 
insulator transition is first-order below Tc ~ 380K [Bj, 
with discontinuities in the ratio of the lattice parameters 
(the two phases being isostructural, however) and in the 
conductivity, accompanied by hysteresis. 

The Mott-Hubbard transition is caused by electron- 
electron repulsion. The fundamental features of this 
transition are traditionally expected ^,^ to be explain- 
able in_tCTms of the half- filled single-band Hubbard 
model 
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which describes electrons hopping on a lattice, interact- 
ing with each other through on-site Coulomb repulsion. 
On the basis of this model the Mott-Hubbard transi- 
tion was studied intensively over the last 35 years. Im- 
portant early results were obtained by Hubbard [g| with 
a Green function decoupling scheme, and by Brinkman 
and Rice |]lO[ with the Gutzwiller variational method m , 
both at r = 0. Hubbard's approach yields a continuous 
splitting of the band into a lower and upper Hubbard 
band, but does not describe quasiparticle features. By 
contrast, the Gutzwiller-Brinkman-Rice approach con- 
centrates on the low-energy behavior, the transition be- 
ing monitored by the disappearance of the quasiparticle 
peak, but does not produce the upper and lower Hub- 
bard bands. A unified approach, treating both the low 
energy and high energy features on equal footing, has re- 



cently become possible within the Dynamical Mean-Field 
Theory (DMFT) |1^, which provides the exact solution 
of the Hubbard model in the limit of infinite dimension- 
ality (or coordination number) |I3]. The complicated 
structure of the self-consistent DMFT-equations makes 
an analytic solution untractable and hence one has to re- 
sort to approximate techniques. In the last few years 
Georges, Kotliar and collaborators performed detailed 
investigations of the metal-insulator transition scenario 
within the DMFT, by employing iterated perturbation 
theory (IPT), exact diagonalization (ED) of small sys- 
tems, quantum Monte-Carlo (QMC) simulations and, at 
T = 0, a projective self-consistent technique (PSCT) 
|0,E3[. While the overall transition scenario reported 
by these authors indeed combines essential features of 
the early approaches, they find the transition to be dis- 
continuous, with hysteresis, due to a coexistence regime 
between the metallic and the insulating phase. For all 
finite temperatures T < Tipt the quasiparticle weight 
disappears abruptly and the gap between the Hubbard 
bands opens discontinuously as a function of U. Hence 
these authors argued that the experimentally observed 
metal-insulator transition in V2O3 can already be un- 
derstood using a purely electronic correlation model. 

At T = 0, numerical renormalization group (NRG) 
studies [Q also found hysteresis and a value for the crit- 
ical interaction of Uc — 5.86 (in our units, see below), 
which agrees with the results of the PSCT ||l^. Never- 
theless, the existence of a preformed gap at T = and the 
corresponding separation of energy scales on which the 
PSCT is based were recently disputed by one of us [[l5[ . 
Finally, a continuous transition with a considerably lower 
Uc was recently observed within the random dispersion 
approximation (RDA (|@). Clearly the Mott-Hubbard 
transition scenario is still very controversial. 

It is the purpose of this paper to carefully re- 
examine the nature of the metal-insulator transition 
within DMFT at finite temperatures. The examination 



of a transition region requires a technique with sufficient 
precision. Since IPT is a rather ad hoc approximation 
scheme, its quahtative and quantitative accuracy is un- 
certain. On the other hand, for T > ED is Umited to 
quite small systems (< 7 sites), so that finite-size effects 
may be considerable. Indeed, although both techniques 
predict the metal-insulator transition to be discontin- 
uous, their quantitative predictions differ substantially 
(see below). In fact, their respective regions of hystere- 
sis do not even overlap much. To resolve these discrep- 
ancies we perform finite-temperature QMC calculations, 
using two different codes to reduce possible systematic 
errors. Although QMC is limited to not too low T- and 
not too large CZ-values, it is still the best understood and 
most thoroughly tested technique presently available for 
the solution of the DMFT-equations. For comparison 
with Ref. ||ll[ we focus on the Hubbard model with a 
semi-elliptical non-interacting density of states (DOS): 
N{E) = y/A- [E/t*Y/{2nt*) if \E\ < 2t* and zero else- 
where. This DOS is realized, e.g., on a Bethe lattice 
with hopping amplitude scaled as i = t*/\/Z, where 
Z — > oo is the coordination number. In the following 
we set t* = 1. In order to study the Mott-Hubbard tran- 
sition, we restrict our calculations to the paramagnetic 
phase and exclude symmetry breaking. The solution of 
the Hubbard model is then determined by the following 
single-site effective action pT| , p"7| -pO[ : 

SoS = - f drdr' J2^Ur)GoHr~r')c,{T') 

In particular, the on-site Green function Gii(T — r') of 
the Hubbard model is identical to the single-site Green 
function G[t—t') = —{Tc {t)c^ {T'))s^in which is imphc- 
itly determined by Scs and the self-consistency relation 
Qi:,{iujri) = [i^n + M — G{i^n)] 1 where Wn = (2n-|-l)7rT. 
The key point in DMFT is thus the accurate calcula- 
tion of the single-site Green function G{t). For this pur- 
pose we use QMC-simulations, which are essentially ex- 
act though computationally expensive ||l^,^ . After dis- 
cretizing the imaginary time into A time slices of length 
Ar — /3/A and performing a Hubbard-Stratonovich 
transformation which introduces auxiliary Ising spins, 
the DMFT-equations are solved by iteration. The num- 
ber of proposed flips of Ising spins ( "sweeps" ) per itera- 
tion will be important here. Each iteration has as input 
the "old" self-energy Y,o\d{i^) and as output a "new" 
self-energy Enow(*w). The rate of change in the itera- 
tion procedure is measured by 77 = A^^ ^^^ \Yio\d{i^n) — 
Sncw(*'^n)l- Experience shows that, for most purposes 
(e.g., calculation of thermodynamic quantities outside 
the critical regime), convergence is reached if 77 < 10~^. 
Physical properties are finally obtained by extrapolation 
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FIG. 1. QMC-resuIts for tlie double occupancy D at 
T — 1/20 using the criterion 77 < 10"'^. The error bars are 
smaller than the data points. All lines are guides to the eye 
only. Inset: Vanishing of the hysteresis as a function of the 
number of further iterations at T = 0.05, f/ = 5 and Ar = 0.3. 



of the simulation results for various Ar- values to the limit 
Ar ^0. 

We found that straightforward implementation of the 
QMC algorithm with 77 < lO^'^ indeed leads to the appar- 
ent convergence of two solutions, i.e., to hysteresis. As 
an example we present results for the double occupancy 
D — (77i|7ii|) in Fig. 1. The upper curve (sohd dots) 
shows results for increasing coupling strength, where we 
use the self-energy calculated for interaction U as in- 
put for the calculation for U + /S.U (here AC/ = 0.3). 
Similarly, the lower curve (open dots) shows results for 
decreasing U . At each [/-value the iteration procedure 
was terminated as soon as the condition rj < 10^'^ was 
satisfied, which typically happens after only a few itera- 
tions. In this way we are able to qualitatively reproduce 
the hysteresis found in Ref. |ll| by means of a standard 
QMC-algorithm with a standard convergence criterion. 
However, as will be shown below, in the case of the 
Mott-Hubbard transition the criterion 77 < 10~^ is not 
sufficient to guarantee convergence to the true solution. 

For this we investigate the stability of the two solutions 
independently under further iterations. As an example 
we consider the solution at T = 1/20 and U = 5. The 
data for Ar = 0.3 are presented in the inset of Fig. 1. In 
order to obtain maximum accuracy we use a large num- 
ber of sweeps per iteration (2 — 3 • 10^). We find that (z) 
both solutions at [/ = 5 in Fig. 1 are unstable, (m) a new 
stable solution is reached after approximately 20 more 
iterations, and (Hi) the solution with the larger number 
of sweeps and time slices (solid dots) reaches equilibrium 
sooner and has smaller fiuctuations. During the itera- 
tion process 77 fiuctuates around (77) ~ 4 • 10^"* and gives 
essentially no information about the distance from equi- 
librium. We obtained similar results for a dense grid of 




FIG. 2. Compressibility k and quasiparticle renormaliza- 
tion factor Z — m/m* vs. U for T = 1/15. All lines are 
guides to the eye only. 



other values for U and At. Extrapolating to Ar -^ and 
combining the results for various C/, we obtain a smooth 
curve D{U) without hysteresis (dashed curve in Fig. 1). 
The same happens at higher temperatures (T = 1/10, 
1/13 and 1/15) and also at lower temperatures, such as 
T = 1/30 and T = 1/35. 

We also studied the quasiparticle renormalization fac- 
tor Z = m/m* and the compressibility n at various tem- 
peratures. After 15-25 further iterations, none of these 
quantities shows hysteresis any more. The results for 
T = 1/15 are shown in Fig. 2. We locate the Mott- 
Hubbard transition at the interaction strength Uc where 
Z{U) and k{U) essentially vanish. The resulting phase 
diagram is plotted in Fig. 3, where the corresponding 
IPT- and ED-result s P J^ and the C/c-values at T = 
obtained by PSCT ppNfRG |l| and RDA |l|] are also 
shown. There is a clear quantitative and qualitative dis- 
crepancy between our numerically exact QMC-data and 
both the ED- and the IPT-results. According to IPT, 
the transition below T — Tipt — 0.088 is discontinu- 
ous. By contrast, we find that the transition from the 
metal to the insulator is continuous down to at least 
T^ 1/35- 0.325 TiPT. 

More information concerning the Mott-Hubbard tran- 
sition can be obtained from the DOS, which we calculated 
from the converged data for the Green function using the 



Maximum Entropy Method (MEM |23|). Our results for 
the DOS as a function of C/ at T = 1/20 are presented 
in Fig. 4. Upon increase of C7, the DOS develops a well- 
defined central peak and shoulders, with the peak pinned 
at its Fermi liquid value. The peak rapidly collapses at 
U ~ Uc (here: Uc — 4.8), leaving an incomplete gap 
p4|. Beyond Uc the system remains semi- metallic un- 
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FIG. 3. Phase diagram of the Hubbard model (paramag- 
netic phase only). Solid squares; continuous metal-insulator 
transition calculated with QMC (error bars include both 
the statistical errors and the uncertainty in Uc due to fi- 
nite temperatures). The dotted line is a guide to the eye 
only. Broad horizontal lines: coexistence region within ED 
[[IIIJ23]. Dashed lines: coexistence region within IPT; the line 
of first-order transitions (full curve) ends at Tjpx (solid cir- 
cle) |11|. The shaded area is a crossover region. Also shown 
are the t/c-values from PSCT/NRG (X) and RDA (0). 



stand the region Uc < U < Ug we studied the tempera- 
ture dependence of the screened local moment (inset to 
Fig. 4) and of the spectrum (not shown) a,t U — 4.8, 
starting from T = 1/20. Upon increase of T the incom- 
plete gap and the screened local moment remain essen- 
tially unchanged. Upon decrease of T the central peak 
and the Fermi-liquid behavior are rapidly restored, and 
the screened moment falls dramatically. The incomplete 
gap and the temperature-independence of the screened 
moment imply that for Uc < U < Ug there are few elec- 
tronic states at the Fermi energy to screen the spins. The 
behavior in this region can be interpreted as emanating 
from a quantum critical point at T = 0, characterized 
by a vanishing DOS at the Fermi level, N{0) = 0. The 
depletion of screening states obtains a natural interpre- 
tation in terms of Nozieres' "exhaustion" scenario p5[ , 
which was recently found to be realized in the periodic 
Anderson model pq |. 

To clarify the influence of band-structure effects we 
also performed calculations for a hypercubic lattice in 
d = 00 including next-nearest-neighbor hopping ampli- 
tudes t' = t'*/^2d{d~l) |^,|8|. Choosing t'* /t* < 
in order to obtain a finite lower band edge, we find the 
hysteresis effects to be strongly suppressed by frustration, 
e.g., no initial hysteresis was observed for t'* jt* < —0.25. 
The phase diagram for i' = 0, previously obtained |2^] 
using QMC and perturbation theory in t* (NCA), is qual- 
itatively similar to that for the Bethe lattice in Fig. 3. 

In summary, we demonstrated that, for temperatures 
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FIG. 4. QMC/MEM-results for the DOS (Bethe lattice) 
for various values of (7 at T = 1/20, and At = 0.017. The 
gap opens at Ug ~ 5.0. The temperature-dependence of the 
screened local moment nis (the product of the temperature 
and the local susceptibility) is shown in the inset for U = 4.8. 



down to Tinin — 1/35 ~ 0.325 Tjpt, the coexistence re- 
gion characteristic of a first-order metal-insulator transi- 
tion disappears in sufficiently careful QMC simulations. 
Further conclusions about the order of the transition can- 
not safely be made. Physically we would expect the tran- 
sition to be a smooth (broadened by finite temperatures) 
but rapid crossover between Uc (where the Fermi liquid 
breaks down) and Ug > Uc (where the gap opens). We 
cannot rule out the existence of a first-order transition at 
even lower temperatures and note that our results may be 
smoothly connected to those from PSCT [Q and NRG 
|l4| at T = 0. However, we found no evidence for the 
existence of a preformed gap at low temperatures, in line 
with the transition scenarios of Refs. ||,|l^,^. - For a 
bandwidth of 0.8 eV pfl] our results imply that down to 



T„ 



70K the metal- insulator transition is continuous. 



Since in the experiment on Cr-doped V2O3 the transition 
is first-order below Tc ~ 380K B, an explanation of the 
nature of this transition apparently requires the inclusion 
of other degrees of freedom, such as lattice [pl|, orbital 
p2| , ^ and possibly higher-spin [|4J effects. 
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